CLNS 97/1516 
CLEO 97-23 



Flavor-Specific Inclusive B Decays to Charm 



CLEO Collaboration 
(October 23, 1997) 

Abstract 



We have measured the branching fractions for B — > DX, B — > DX, and 
1? — ► DXl + v, where is an average over 5° and "D" is a sum over 
D° and -D + , and "D" is a sum over Z)° and D~. From these results and some 
previously measured branching fractions, we obtain B(b — > ccs) = (21.9 ± 



3.7)%, B(b -» as) < 6.8% @ 90% c.l, and B(D° -> if-7T+) = (3.69 ± 0.20)% 



0. 



Implications for the "S semileptonic decay problem" (measured branching 
fraction being below theoretical expectations) are discussed. The increase in 
the value of B(b — > ccs) due to S — > eliminates 40% of the discrepancy. 

PACS numbers: 13.25Hw, 14.40.Nd 
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I. INTRODUCTION 



There has been a longstanding problem in heavy flavor physics of the measured B semilep- 
tonic decay branching fraction |l] being smaller than theoretical expectations [0,0. One 
possible explanation is a larger-than-expected flavor-changing neutral current (FCNC) 
contribution, due to new physics. Another || is an enhanced rate for b — > ccs' (s' denotes 
the weak isospin partner of c). An argument against an enhanced b —>■ ccs' rate is that it 
would conflict with the measured branching fraction for B — > DX plus B — > DX. That 
measurement relies on a knowledge of B(D° — > K~tt + ), however, and if that is in error, the 
measurement of the branching fraction of B to charm or anticharm will be also. We address 
all three issues by measuring the yields of the flavor-specific inclusive B decay processes 
B — > DX, B — > DX, and B — > DX£ + u in a sample of BB events in which at least one B 
decays semileptonically. (Herein, "5" represents an average over B° and B + , U D" a sum 
over D° and D + , and "D" a sum over D° and D~ [|J. We use the term "upper vertex D" for 
a D produced from the charm quark from W — > cs, and "lower vertex D" for a D produced 
from the charm quark from b — > c.) 

These yields, and ratios among them, provide information on the above-mentioned issues 
as follows: 

(i) The fraction of semileptonic B decays that proceed through B — > DXi + v, f$L, 
differs from 100% only because of small contributions from b — *> uiv and B — > D~KX£ + i> 
("lower vertex D s "). The measured fraction is inversely proportional to the assumed D 
absolute branching fraction (in our case B(D° — > K~7r + )) and scaling the yield to agree with 
expectations gives a new method to measure that branching fraction. 

(ii) The fraction of all B decays that proceed through B — > DX, f a u, differs from 100% 
because of b — > u decays, lower vertex D s , formation of cc bound states, formation of charmed 
baryons, and FCNC processes such as b — > sg, b — > dg, b —>■ sqq, b — * dqq (which we will 
refer to collectively as "6 — > sg"). As all processes except b — > sg have been measured, the 
ratio faii/fsL provides a measurement of the branching fraction for b — > sg. By taking the 
ratio of f a u to fsi, rather than just using f a u, we eliminate dependence on the D° — > K~n + 
branching ratio, and reduce dependence on D detection efficiency. 

(iii) The process B — > DX proceeds via the the quark-level process b — > ccs', and thus 
the ratio of the yields for B — > DX and B — > DX, i.e., ratio of upper to lower vertex charm, 
provides information on the rate of that process relative to b — > cud'. 

The typical inclusive B decay branching fraction measurement averages over B and B 
initial states for a given final state, and consequently averages over particle and antiparticle 
final states for a given initial state (B or B), losing the flavor-specific information sought here. 
In 1987 CLEO developed a technique for measuring inclusive B decay branching fractions 
separately to particle and antiparticle final states, and applied it to inclusive kaon decays 
|5|||. Here we apply similar techniques to inclusive charm decays. 

The principle underlying the 1987 technique is that if one B from a BB pair from the 
T(4S) decays semileptonically, with a high momentum lepton, then the other decay products 
from that B will have substantial angular correlations with the lepton, tending to come 
off back-to-back to it, while the decay products from the other B have negligible angular 
correlations with the lepton. The lepton tags the flavor of its parent B, and thus also the 
other B (with a correction needed for mixing). By plotting the distribution in the angle 
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between D£ + (and D£~) pairs, and separately the distribution in the angle between D£~ 
(and D£ + ) pairs, and extracting an isotropic component and a peaking component from 
each, yields are obtained for four processes: B — > DX£ + v, B — > DX£ + v, B — > Z)X, and 
5 -> DX. Of these, B -> Dll+y should be zero. 

For low Z) momenta, the technique just described loses statistical power and becomes 
sensitive to the shape assumed in fitting for the peaking component. (In the limit that the 
D momentum vanishes, the D-lepton angular correlation clearly contains no information.) 
Consequently, we have developed a second technique, based on charge correlations alone. 
We measure three yields: the number of D£~ (and D£ + ) pairs, equal to the sum of B — > 
DXl + v and B — > DX yields in the lepton-tagged data sample; the number of D£ + (and 
D£~) pairs, equal to the sum of B — ► DX£ + v and B — > DX yields in the lepton-tagged 
sample; and the number of D (and D) mesons in an untagged sample, equal to the sum of 
B — > DX and B — > DX yields in the untagged sample. Using the fact that the rate for 
B — > DX£ + v vanishes, and scaling the last-mentioned yield by the ratio of the sizes of the 
tagged and untagged data samples, these yields give the yields for the other three processes: 
B — > DX£ + u, B — > DX, and B — > DX. Using a combination of the angular correlation and 
charge correlation techniques, we have obtained those three yields for the sum of D° and D + 
mesons. 



II. PROCEDURES 

The data were taken with the CLEO detector |7j at the Cornell Electron Storage Ring 
(CESR), and consist of 3.2 fb _1 on the T(4S) resonance and 1.6 fb~ x at a center-of-mass 
energy 60 MeV below the resonance. The on-resonance sample contains 3.3 million BB 
events and 10 million continuum events. The CLEO detector measures charged particles 
over 95% of 47r steradians with a system of cylindrical drift chambers. Its barrel and endcap 
Csl electromagnetic calorimeters cover 98% of 4tt. Hadron identification is provided by 
specific ionization (dE/dx) measurements in the outermost drift chamber and by time-of- 
flight counters (ToF). Muons are identified by their ability to penetrate iron; electrons by 
dE/dx, comparison of track momentum with calorimeter cluster energy, and track/cluster 
position matching. 

We select hadronic events containing at least 4 charged tracks. We require a value 
of the ratio of Fox- Wolfram parameters | ^ e H 2 /H < 0.5, to suppress continuum 
events. Events containing at least one lepton with momentum between 1.5 and 2.8 GeV/c 
and surviving a ip — > £ + £~ veto are scanned for D°, D + , and charge conjugates. (For the 
untagged sample, we drop the lepton requirement.) We detect D° and D + via the K~n + 
and K~"-n + it + decay mode, respectively. Tracks used as candidate D decay products must 
have dE/dx and/or ToF values within 2a of expectations for the particle assignment made 
(K or 7r). For D° — > K~tt + , particle identification must rule out the D° — > tt~K + option. 

For candidate Z)'s, we histogram the Kir (Kmr) mass for four intervals in cosOn-e and 
four intervals in D momentum, separately for the two charge correlations with the lepton. 
These 64 mass distributions are fit to double- Gaussian signal peaks and polynomial back- 
grounds, to extract D yields. These are corrected for detection efficiency, determined by 
Monte Carlo simulation augmented by studies of particle ID efficiency that use data (a sam- 
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pie of D* + — > D°7T + , D° — > K~n + events). We perform small subtractions for continuum 
background (using below- Y(4S)-resonance data) and for hadrons misidentified as leptons 
(using hadrons in place of leptons and weighting by the "faking probability"). Small correc- 
tions are made to the D° yields for the singly- Cabibbo-suppressed decays D°/D° — > K~K + 
and D°/D° — > tt~tt + which combine with a single failure of particle ID to make satellite 
peaks, for the doubly-Cabibbo-suppressed decay D° — > K + n~ ||, and for double failures of 
particle ID, with tt~K + treated as K~7t + . A small correction is made to D + yields for the 
decay — > K~K + tt + with the K + misidentified as a 7r + . 

The D yields for each momentum interval, charge correlation, and D type, are his- 
togrammed vs. cos#d_^, 16 distributions in all. For the high D momentum intervals 1.3 
- 1.95 and 1.95 -2.6 GeV/c, we fit the £~D angular distributions to an isotropic compo- 
nent and a backward-peaking component, with fitting functions obtained from Monte Carlo 
simulation. We fit the £ + D angular distributions to an isotropic component alone. For the 
low D momentum intervals 0.0 - 0.65 and 0.65 - 1.3 GeV/c, we use the charge correlation 
technique, summing over cos8d-i- We sum the yields so obtained over D momentum in- 
tervals, and over charged and neutral ZTs, correcting for D° and D branching fractions, 
using B(D° -> K~TT + ) = 3.91% |LT|, and B(D+ -> K~TT + Tr + )/B(D° -> K-n+) = 2.35 
||1 1|| . We obtain yields for D and lepton from the same B, and from different i?'s, as follows. 
N(D£~+D£+, sameS) = (3. 75±0. 11) x 10 5 , N(D£~+D£ + , different B's)= (6.66±0.77) x 10 4 , 
and N(D£+ + D£~, different B's)= (3.18 ± 0.08) x 10 5 , in a sample containing 4.24 x 10 5 
leptons. For illustrative purposes, we show cos8d_£ distributions summed over momentum 
intervals and over D° and D + , as Fig. 1. The £~ D + £ + D distribution shows strong back- 
to-back peaking from B — > DX£ + u, while the £~D + £ + D shows no such peaking, due to 
the nonexistence of B — > DX£ + v. One also notes a much larger isotropic component in 
£r D + £ + D, because of the large rate for B — > DX and a small rate for B — > DX (and a 
small rate for mixing B° — > B° — > DX). 

If the lepton and D come from the same B, then the lepton tags that B correctly. The 
lepton can't be from decay of D, because the D was detected via a hadronic decay mode. 
It can't be from if), because the rate for B — ► ipDX is negligible. If there are two D's from 
the same B, leptons from either will be below our 1.5 GeV/c momentum cut. If the B has 
mixed, nonetheless the lepton correctly tags the b flavor at the instant of decay, which is 
what is relevant for understanding the D from the same B. But, if the lepton and D come 
from different i?'s, then the tagging of both i?'s is now imperfect: the ancestor of the lepton 
because leptons from charm decay and leptons from ip now contribute; and the ancestor of 
the D for those reasons and in addition because of B° — B° mixing. Corrections are thus 
required when using the yields involving lepton and D from different -B's. These corrections 
depend on f m (the probability that a lepton mistags its ancestor B), and x (the mixing 
parameter) . 



III. RESULTS AND INTERPRETATION 

We extract three distinct pieces of physics from the three yields given above. For each, we 
have considered systematic errors due to uncertainties in each of the previously-mentioned 
corrections, uncertainties from fitting mass peaks and cosOd-i distributions, and uncertain- 
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ties in efficiency and D branching fractions. 

(i) First consider T(B — > DX)/T(B — > DX), the ratio of "upper vertex" charm to "lower 
vertex" charm. This ratio U / L is obtained from x = N(D£~ + D£ + , different S's) /N(D£ + + 
D£~ , different S's) , by correcting for mixing and mistags. U/L = (x — F m ) / ( 1 — xF m ) , where 
F m = (fm + f')l (2 - f m - /'), and f = f m + X~ V™X- We use x = 0.157 as measured 
by CLEO with dileptons [|12|], and f m = 0.027 as found there, thereby achieving cancellation 



of some systematic errors in F m , giving F m = 0.112±0.011. From the yields given above, 
x = 0.210 ± 0.025, leading to 

S| 2 Dx] = °- 100 ± °- 026 ± °- 016, (1) 

where the first error is statistical and the second is systematic, dominated by the uncertainties 
in mixing correction(± 0.012) and the cos#d_£ fitting function (± 0.008). This result is 
surprisingly large, as conventional wisdom held that b — > ccs would hadronize dominantly 
into D s . However, Buchalla et al. @ have argued that the D°, D + component should be 
substantial. 

In Fig. 2 we plot the momentum distribution of these "upper vertex" D°, D + , obtained 
by applying the analysis just described to each of the four D momentum bins. The spec- 
trum is softer than that for "lower vertex" Z)'s, also shown. It is well described by 3-body 
[)(*)£){*) Ki*) phase space, if one allows one or two of the particles to be the vector states. 
CLEO has observed such decay modes |T3| . 



(ii) Next consider the fraction of all B decays to D, f a u, divided by the fraction of semilep- 
tonic B decays to D, fsL, i-e., the double-ratio of widths / ^B^xe+J) • We obtain 

this from the ratio of yields N(D£ + + D£~, different B , s)/N(D£~ + D£ + ,s&meB) = z raw . 
Corrections are required to the "different -B's" yield for mixing and mistags. Also, lep- 
tons from unvetoed ip and from secondary decays (3.3 ± 0.7% of all leptons) don't con- 
tribute to the peaking yield, and so a correction is required for that, leading to z cor = 
0.967z raw /[(l - 0.5/ m - 0.5/')(l + F m U/L)}, where U/L = 0.100, as found above. Applying 
all corrections, we have 

T ^l T ~^m-^ -,0 1± 0.03 4± 0. 01 , (2 ) 

One expects both f a u and fsi to be close to 1.0. The first ratio will be less than 1.0 because 
of b — » u transitions (2| V u b/V c b\ 2 , where the 2 is a phase space factor), lower vertex D s (2%), 
bound cc states (3.0 ± 0.5% 10), baryons (6.5 ± 1.5% [fTSfl), and b — >• sg (to be extracted). 
The second ratio will be less than 1.0 because of b — > u transitions (3| V u b/V c b\ 2 , enhanced by 
the 1.5 GeV/c lepton momentum requirement), and lower vertex D s (1.0 ± 0.5%, suppressed 
by the lepton momentum requirement). These lead to 

faii/fsL = 1-0 + | V ub /V cb \ 2 - (0.010 ± 0.005) - (0.030 ± 0.005) - (0.065 ± 0.015) - B(b -> sg) 

Here b — >• sg is symbolic for all FCNC processes. Using \ V u b/V c b\ 2 = 0.008 ±0.003, we obtain 
B(b — > sg) = (0.2±3.4±1.5±1.7)%, where the first error is statistical, the second systematic 
on z, and the third the uncertainties in Expression (3). From this we obtain an upper limit 
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B(b — ► sg) < 6.8%, @ 90% CL. The dominant components of the systematic error on z are 
from mixing (±1.2%) and unvetoed and secondary leptons (± 0.6%). 

(hi) Finally consider the fraction of semileptonic B decays to D° or D~, i.e., fsL = T(B — > 
DXl + v) /T(B -> X£+z/). We obtain this fraction by dividing the yield N(D£-+D£ + , sameS) 
by the number of leptons from B semileptonic decay, 96.7% of the total of 4.24 x 10 5 leptons 
in our sample. We find 0.914 ± 0.027 ± 0.042. This number is inversely proportional to 
the value used for B(D° — > K~tt + ). The expected value of the ratio of widths is T(B — »■ 
DX£+u)/T(B -> X£ + v) = 1.0 - 3\V ub /V ch \ 2 - 0.010 ± 0.005(forfi -> D+KX£~v). Taking 
3|WKfc| 2 = 0.023 ± 0.008, we find the expected ratio of widths to be 0.968 ± 0.010, 
differing from the measured value by one standard deviation. We set measured and expected 
values of the ratio equal to each other, and solve for the D° branching fraction, finding 
B(D° -> K-tt + ) = (3.69±0.11±0.16±0.04)%, where the first error is statistical, the second 
systematic in the measured ratio, and the third systematic in the predicted ratio. The 
dominant systematic errors are from uncertainties in D detection efficiency (± 0.10%), mass 
peak fitting (± 0.09%), and the ratio of D + to D° branching ratios (± 0.08%). This value 
for the branching fraction, (3.69 ± 0.20)%, is to be compared with recent measurements by 
CLEO of (3.91 ± 0.19)% || and (3.81 ± 0.22)% |T§, by ALEPH of (3.90 ± 0.15)% ]T7f 
and the PDG value of (3.83 ± 0.12)% |J[. 

IV. THE B SEMILEPTONIC DECAY BRANCHING FRACTION PROBLEM 

In Table I we list all the components of B decay, give their branching fractions (based on 
measurement or theory), and see if they sum to 100%. We express some in terms of bsL, the 
B semileptonic decay branching fraction, for which we use |1] (10.49 ± 0.46)%. The factor 
of 0.25 for b — * (c or u)rv is a phase space factor. The factor r u & for b — > (c or u)ud' would 
be 3 from color counting, but with QCD corrections |19| is 4.0 ± 0.4. This analysis has two 
pieces of information to add to Table I. First, the upper vertex D°, D~ contribution of (7.9 
± 2.2)% is obtained from our measured value of T(B D° or D + X)/T(B D°otD-X) 
combined with the rate for inclusive D° + D + (63.6% + 23.5%) |[20|| , and leads to a branching 
fraction for b — > (coiu)cs' of (21.9 ± 3.7)%. Second, we have a value (with large errors) for 
the FCNC term. One sees that the upper vertex D°, D~ contribution accounts for close to 
half of the shortfall of the sum of all modes from unity. The remaining shortfall is less than 
2 standard deviations. If we adjust r u d to bring the sum to 100%, we find r u d = 5.2 ± 0.6. 
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FIG. 1. Yield of D£ events vs cos 9 D - £ . D°e~ + D+£- plus charge conjugate, summed over 
D momentum are shown as solid circles, while D°£~ + D~l~ plus charge conjugate, summed over 
D momentum are shown as open squares. 
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FIG. 2. D momentum distributions. "Upper vertex" D° + D + , i.e., from B — > DX are shown as 
solid squares, while "lower vertex" D° + D + , from B — > are shown as open squares and "lower 
vertex" D° + D + , from B — > DXiu are shown as solid circles. Vertical scale gives branching fraction 
per unit momentum, for upper and lower vertex Z)'s, and same divided by total semileptonic decay 
branching fraction for semileptonic D's. 
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TABLE I. All components of B decay, with their branching fractions. Upper vertex D° and 
D~ , and b — > s/d g, s/d qq, are from this analysis. The branching fractions for the separate 
components making up b — > (c or u) cs' are shown parenthetically. 



b decay modes 




Branching fraction (%) 




b — ► (c or u) 






10.5 ± 0.5 


6 — > (c or u) iiv 


^SL 




10.5 ± 0.5 


6 — > (c or n) Tf 


0.25&5L 




2.6 ± 0.1 


b — > (c or u) wc?' 






42.0 ± 2.0 ± 4.2 


6 — > (c or n) cs' 






21.9 ± 3.7 






(10.0 ± 2.7) 




(cc) 




(3.0 ± 0.5) 




baryons 




(1.0 ± 0.6) 




upper vertex D°, D~ 




(7.9 ± 2.2) 




b — ► s/d <?, s/d qq 






0.2 ± 4.1 


TOTAL 






87.7 ± 7.4 
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